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Abstract 



i-i^ . We examine SU{2) gauge theory in 3 + 1 dimensions at finite tem- 

perature in the vicinity of critical point. For various lattice sizes in time 

#vj , direction (A'^^ — 1,2,4,8) we extract high precision values of the inverse 

^ ' critical coupling and critical values of the 4-th order cumulant of Polyakov 

00 ' loops (Binder cumulant). We check the universality class of the theory 

^^ , by comparing the cumulant values to that of the 3D Ising model and find 

C^^ ' very good agreement. 

^^ I The Polyakov loop correlators for the indicated lattices are also mea- 

sured and the string tension values extracted. The high precision values 
of critical coupling and string tension allow us to study the scaling of di- 

t*^ , mensionless T^j ^fa ratio. The violation of scaling by < 10% is observed 

c I ^ ■ as the coupling is varied from weak to strong coupling regime. 

>: 

^ ■ 1 Introduction 



The original motivation for this study was to re-derive the values of critical 
couphng for finite temperature SU{2) gauge theory formulated on lattices with 
different time-like extent Nr at the highest precision possible with modern com- 
putational resources (we utihze a small 20-30 node computer cluster). For this 
we rely on the property of universahty of SU{2) gauge theory at the second 
order critical point. We adopt a standard for such models procedure of locat- 
ing the position of critical point, which amounts to measuring the 4th order 
(Binder) cumulant Ij 54 for Polyakov loops P 

54 = 1-^, ^=]^E^Trnf/.,,,o (1) 



on Nt X A^^ lattices in the vicinity of phase transition. Note that we use the 
original form of the Binder cumulant which differs by a constant factor 1/3 from 
the normalized version frequently used in lattice gauge theory literature. The 
finite size scaling (FSS) of 34 with lattice size N^ in the vicinity of the critical 
point is knownj^l [5] 

9i.N„ ~ .94,00(1 + a^tNl''' + a2N-y' + • • •) (2) 

where t = {T — Tc^oq)/Tc,oo is the reduced temperature and j/i = —lo > is the 
exponent of the largest irrelevant scaling field. 

The scaling of the temperature value of the intersection point of Binder 
cumulant curves 34 (i) on A^^r and N'^ = bN^, b > 1 lattices is 

t* = _^iv-^.-i/-i:i^. (3) 

ai " 1 - foi/^ '^ '' 

It is convenient to define Lb = Na^'~'^^''{l-b-y')/{l-b^/''), so that jS]) becomes 
a simple linear function of Lj,. Substituting t* into ^ results in the scaling for 
the intersection 174 value 

9t = 54,00 [l + a,N-y^ x_^-il. ) ■ (4) 

Here we define Lg = iV-^i(l - 6~^i"i/'')/(l- 6"^/''), so that Q is linear in L^. 
The well-known renormalization group relationship allows us to connect the 
lattice spacing to the coupling /3 = 2Nc/g'^, [Nc = 2) 

Using the fact that the temperature is T = l/{Nra) one obtains the following 
expression for the reduced temperature in the vicinity of the transition point 



In the studies of Nr — 8 and 16 phase transitions by Fingberg et al.[3] no 
deviations from the leading term were found. Therefore one can assume linear 
correspondence between the reduced temperature and reduced inverse coupling 
and perform the scaling studies in terms of the inverse lattice coupling /3. 

In their renowned work Svetitsky and Yaffe[5] conjectured that the univer- 
sality class of the d + 1-dimensional SU{2) gauge theory is the d-dimensional 
Ising model. Since then it was confirmed in numerous numerical simulations, 
e.g. Ref. [5] for N^- — 4 theory. Indeed, we observe that for all lattices in 
the thermodynamic limit the g^ curves intersect in the vicinity of the 31? Ising 
value, albeit with different degree of accuracy. 

In section[2]we measure values of Binder cumulant in the vicinity of the crit- 
ical point and determine critical values of inverse coupling and Binder cumulant 



for Nr — 1, 2, 4, 8 lattices. Our treatment of different Nr lattices is not uniform. 
We conduct the FSS study for iV,- = 4 and 2 lattices, while relatively large N„ 
values used for Nr = I studies prevented us from observing a significant scal- 
ing behavior. The Nr = 8 lattice allows for the FSS study, however it is very 
expensive to simulate and therefore is studied assuming the Ising universality 
class. 

The small Nr = 1 and 2 lattices are of special interest in decimation studies, 
since the iterative block-spinning procedure has to be stopped when the smallest 
(or next to smallest) lattice is reached. Therefore for N^ = 1 lattice we consider 
a specific lattice formulation suitable for decimation. 

In section [3] we study the quark-antiquark static potential and extract the 
string tension for Nr/i > 2 lattices at temperature T — Tc/3. The knowledge 
of critical couplings allows us to construct the dimensionless ratio Td \fo and 
study its scaling with coupling. We observe relatively small < 10% violation of 
scaling at strong coupling /3 ~ 1.8738 value. 

2 Finite temperature phase transition: Monte 
Carlo study 

In all simulations performed in this study we use the standard Wilson action. 
Per one updating sweep we perform two overrelaxation steps and one heat- 
bath update. We use a standard acceptance improved heat bath updating 
procedure[7l [8]. Measurements are performed every 2 — 20 sweep. We use 
10 — 80 independent runs (at different initial random generator seeds), each run 
is averaged into a single bin. This allows us to gain better statistics and avoid au- 
tocorrelations. All errors are computed with the jack-knife method with respect 
to these bins, except when indicated differently. For various coupling values 
/3 = 4/(/^ in the visinity of the finite temperature confinement-deconfinement 
phase transition after initial equilibration we measure Polyakov loops and com- 
pute the Binder cumulant. The number of sweeps needed for the system to 
reach equilibrium is estimated by observing the Monte Carlo time evolution 
of the Polyakov loop and plaquette estimates for each of the lattice sizes and 
typical equilibration times are > lO'^ sweeps. 

2.1 A^^ = 4 lattice 

We start with perhaps the most studied lattice TV^ = 4. At this point we 
do not need to assume any particular universality class. We simulate Na = 
8, 10, 16, 20, 24 and 32 lattices, with typical statistics of 30 x 40000 configura- 
tions. 

The results for the Binder cumulant measurement are presented in Fig. [T] 
together with the fitting lines .g4(/3) = a/3 + 6 (see Tab. [l]for fitting parameters). 
The data is fitted to a straight linqj in the vicinity of the transition point with 
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Figure 1: The Binder cumulant 34 for iV^ = 4, 7V^ = 8,10,16,24,32 lattices. 
Linear fits are represented by solid lines, the T^^ estimate with errors by vertical 
hnes. 



a standard x^ procedure. The errors in the fitting parameters a and 5 are 
obtained by error propagation. The resulting goodness of fit in all cases is 
Q « 0.6 — 0.9, which confirms that at the considered /3 values the 54 curves are 
well approximated by lines. In Tab. Owe list the values /?* and g\ of intersection 
points for different pairs of hnes (TVo-, N'^ = bN^)- 

First we check the assumed scaling behavior by plotting the rescaled Binder 
cumulant g^^Na {(3— (3c) / Pc) for various lattice sizes, see Fig. [2 For the critical 
coupling value we use /3c = 2.2991, which we will obtain later in the FSS study. 
Indeed, we observe that the curves fall on top of each other. The only noticeable 
deviation from the scaling behavior can be observed for the smallest considered 
lattice Na = 8 far away from the transition point. This is obviously due to the 
effects from next to leading order scaling terms. The important observation is 
that in the region where we perform linear fits there is no deviation from the 
scaling for all the lattices and also apparently the linearity holds. 

Next we study the scaling of the pair-wise intersection point coordinates 
according to ([3]) and Q. First we look at the (3* coordinate of intersection 
points versus L^, see Fig. [31 and then at the gX coordinte of intersection points 
versus Lg, see Fig. 21 For this we need to know the values of critical exponents 



^ We have found that the reweighting of Polyakov loops to new /3 values has extremely 
short range. Therefore we do not use reweighting and rather use the linear fitting. 



Na_ 

a 
b 



10 



16 



20 



24 



32 



4.17(29) 5.77(36) 11.80(56) 17.6(1.2) 20.6(1.3) 34.7(1.9) 
-9.13(66) -12.89(82) -26.7(1.3) -39.9(2.8) -46.9(2.0) -79.3(4.3) 



Table 1: Resulting parameters of the x^ fit of the Binder cumulant to a linear 
gi{(3) = a(3 + b function. 
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Figure 2: The rescaled Binder cumulant gi{Na {(3 — (ic)/ fie) for Nr = 4, Na 
8, 10, 16, 24, 32 lattices. 
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u and j/i . As we will see (from the quality of data) it is not practical to extract 
them from the data, instead we take them to be equal to the Ising values known 
with good accuracy 1/v — 1.5887(85) and yi — 0.812 [9]. For a self-consistency 
check we verify that the scaling resulting from use of these values is adequate. 
For each smaller lattice size from N^ = 10, 16, 20, 24 set (represented by 
a different point type/color on the figure) we take various possible b values, 
which index the larger lattice. As one can see from Fig. [3] despite the fact that 
the simulation statistics is suflricient for the precise location of the intersection 
point the error bars can become large as intersecting lines approach the coUinear 
limit. We performed a linear fit in accordance with ^. The goodness of fit 
is Q = 0.66 which suggests that the scaling equation is correct. The limit 
Lb —^ Q corresponds to the thermodynamic limit and yields /3c, oo = 2.2991(2). 
Note that the largest lattices intersection (Nr = 24 and N^ = 32) happens at 
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Table 2: Pairwise intersection coordinates in g4 — (3 plane: the inverse lattice 
coupling (top) and the 4th cumulant (bottom); N^ — 4. 
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(3* = 2.2993(3) which agrees with the thermodynamic limit result. These results 
should be compared with the value Pc = 2.2985(6) of Ref. ^ or 2.2986(6) 
(intersection 12 and 18) of Ref. [TU] and indeed we find a very good agreement. 

The value of Binder cumulant for the 3D Ising model is well known. In the 
lattice gauge theory literature the early estimate of Ferrenberg et al. 54"*"^ = 
0.470(5J3 [9] is often used. In this work we are using the more recent and 
more accurate estimate of Hasenbusch et al. [IT Q =< m^ >^ / < m^ > = 
0.62393(13) (here m is the Ising model magnetization), which translates into 
0.46575(11). 

The data of Tab. [2] and Fig. [1] strongly support the fact that the Binder 
cumulant reaches the Ising value for the larger lattice intersections. Therefore 
it is safe to assume that within the accuracy of this study the Nt = 4 model is 
indeed in the 3D Ising model universality class. 

We continue the study of the scaling of pair-wise intersection point coordi- 
nates with the size of participating lattices. Now similarly to /3* analysis we 
plot the gl values of intersection points versus Lg , see Fig. |4l As one can see 
the data does not allow to perform a quality fit especially in the thermodynamic 
limit Lg -^ 0. Since we already assumed that the Ising universality class holds 
we may fit 34 directly to a* f] + 54'*™^. The reduction of the degrees of freedom 
allows for a better defined fit. The goodness of fit is Q = 0.80 meaning that the 
scaling function is plausible. Also it represents a self-consistency check for the 
assumption of the Ising universality class. 

At this point it is worth to mention that in the studies of the Dyson's hi- 
erarchical model [H] it was observed that a too coarse scale in the /? sampling 
may lead to inaccurate locations of the intersections when linear fits are used. 
This inaccuracy may generate sizable errors in the Lj, —f and Lg —f limits. 

As a side note we would like to point that insufficient statistics prevented 



^The error is from the published figure. 
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Figure 3: The /3* coordinates of the pair- wise intersection points for N^ = L ^ 
10, 16,20,24 and various h lattices. The sohd hne represents a hnear fit for all 
the presented data. Nr — 4. 



us from performing the fits inside separate bins and thus obtaining the relevant 
errors by use of the jack-knife method on this bins. Instead the fits were per- 
formed on a whole set and the errors were obtained by simple error propagation, 
which may undermine their accuracy. 



2.2 Nr = 2 lattice 

Next we study Nt = 2 lattice. This lattice is somewhat special because unlike 
larger Nr lattices as a result of periodic boundary conditions in T direction 
every pair of spatial links from different time slices is connected only by a single 
time-like link. 

For this lattice we generated on average from 10 x 0.5 * 10^ configurations 
for No- = 16 to 10x0.3*10^ configurations for iVo- =32 lattices. The results for 
Binder cumulant measurements in the vicinity of transition point are presented 
in Fig. [5l As one can see from the figure the linear fit lines intersect clearly 
below the Ising value and at what appears to be a single point. We collect 
the coordinates of the intersection points of the fitting lines in Tab. [31 We 
observe that for the intersection coordinate f3* of larger lattices there is very 
insignificant scaling change. Therefore for this quantity we do not perform the 
FSS analysis and take as the transition coupling the value of intersection of lines 







g 

Figure 4: Same as Fig. [31 but for the Binder cumulant. The horizontal hne is 
the Ising 54^™^ = 0.46575 value. 



of the largest 32 and 24 lattices: /3c = 1.87348(2). 

For the other intersection coordinate gl there is a noticeable scaling behavior, 
therefore we perform the FSS analysis analogous to the presented earlier {Nr = 
4). The results are plotted in Fig. [S) It is interesting that the smaller N^ = 10 
lattice is not consistent with the scaling behavior, while it was for A^^ — 4 
lattice. This can indicate that the sub-leading effects can be stronger here. 
For the larger lattices, however, the fit is good and indicates that the assumed 
3D Ising universality class is correct. Also the assumption of the 31? Ising 
universality class is suported by independent studies of critical exponent v |13| . 

It is expected that for lattices with smaller Nt it is sufficient to consider 
lattices with correspondingly smaller No- values, therefore it is surprising that 
the intersection point of No- — 32 and 24 lattices is located statistically signifi- 
cantly below the Ising g^^"''^ value and one needs even larger lattices to reach 
the value corresponding to the thermodynamic limit. 

For Nt = 2 we found several estimates in the most recent literature all in 
good agreement with our result. On 2x 12^ lattice the estimate is /3c = 1.877 (no 
error given, presumably 1.877(1)) [14], while more recent estimate from 2 x 30^ 
and 2 X 40^ lattices gives /3c = 1.8735(4) [HITS]. 
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Table 3: Pairwise intersection coordinates in g^ — (3 plane: the inverse lattice 
coupling (top) and the 4th cumulant (bottom); N^ = 2. 

2.3 A^^ = 1 lattice 

Here we look at even more exceptional Nr — 1 lattice. This lattice is very special 
and allows for two formulations. One formulation is motivated by considering 
the sum over plaquettes in the partition function. It is natural to assume that 
there is one time-like plaquette for each space coordinate. This case requires 
special treatment of space-like links since they have conjugate staples only in 
one time direction (we choose positive direction) unlike two directions (positive 
and negative) for other directions. The other formulation can be motivated 
by considering the Nj- = 2 lattice and performing factor 2 decimation in time 
direction (removing one time slice of links). In this formulation there is no 
special treatment of conjugate staples of space-like links. One has to consider 
staples in both positive and negative direction for each direction (including the 
time direction). 

We start with the first formulation of N^- = 1 lattice, which we call time-like 
plaquette single counting. Afterwards we will consider the second formulation 
which we refer to as time-like plaquette double counting. 

For single counting formulation we present the results in Fig. [T] Similarly 
to Nr — 2 case we do not observe any scaling behavior since the lines cross at 
a single point. The linear fitting lines intersection coordinates are: 
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= 0.85969(13),. g^ = 0.4535(20) 
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= 0.85989(11), g^ = 0.4573(28) 
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= 0.85997(10), g^ = 0.4606(51) 



Note little change in the /3* and g| coordinates, there is virtually no scaling 
to be extracted here. The intersection of the largest lattices (32 and 24) defines 
the critical coupling /3c = 0.85997(10). The Binder cumulant is within one sigma 
from the Ising value. 

Next we study the second formulation with time- like plaquettes double count- 
ing. We summarize the results in Fig. [51 The pairwise fitting lines intersection 
coordinates in this case are: 
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Figure 5: The Binder cumulant 94 for Nr ^ 2, N„ = 10,16,24,32 lattices. 
Linear fits are represented by solid lines, the Tc estimate with errors by vertical 
lines. 



N^ = 16, 7V; = 24 - 13* ^ 0.86228(6), gl = 0.4606(12), 
7V^ = 24, iV; = 32 - /3* =0.86226(6),. g4 = 0.4598(31). 

Again there is almost no scaling behavior in the coordinates. The values of 
Binder cumulant arc comparable to the previous formulation, however the in- 
verse critical coupling is shifted toward larger values. Also we observe that the 
Binder cumulant deviation from the Ising value is small (< 2a away). 

The only known result in the recent literature for A^,- — 1 lattice is /3c = 
0.8730(2) [T7]. However, the authors do not indicate how it was obtained. 



2.4 A^^ = 8 lattice 

Next we look at A^,- = 8 and A^^ = 16, 24 and 32 lattices. We present the results 
in Fig. [9l For lattice N^ = 16 we performed 10^ sweeps, while for No- = 24 
and 32 lattices we performed 0.2 x 10^ sweeps (measuring every 20) at the three 
beta values closest to the transition. 

The approach adopted for smaller lattices is not very efficient here. The 
fitting lines for larger two lattices intersect at 
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/3 = 2.5113(19), 34 = 0.472(13). 
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Figure 6: Same as Fig. [H but for Nr = 2 lattice. 



As one can see the errors are quite significant here and better statistics is needed. 

On the other hand if one assumes the universalhty class of the 3-D Ising 
model it is possible to obtain better results. Our strategy is similar to one used 
in Ref. [TH] ■ We look at the points of intersection of linear fitting curves with 
the 3-D Ising Binder cumulant value 

iV„ = 16 f3^ 2.5143(0) 
iV^ = 24 (3 = 2.5104(2) 
iV^ = 32 /3 = 2.5107(8). 

The crossing of the fitting line with the Ising value defines the transition 
point. Note that N^ = 16 lattice intersects with N^r = 24 slightly off, while N^r = 
24 and N^ = 32 at the Ising value. This indicates that the thermodynamic limit 
has set in for N^ > 24. The uncertainties of the fitting parameters contribute to 
the error of the critical coupling. Possibly the best result is for N^ = 24 lattice 
since the larger lattice is a bit noisier, therefore we define (3c = 2.5104(2). 

We compare our result to a similarly performed study of Ref. [T5], which 
is /3c = 2.5105(10) for iV^ = 8 and N^- = 32,40,48 lattices. The Gaussian 
difference test yields Q = 0.92, indicating good correspondence of the results. 
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Figure 7: The Binder cumulant g^ for A^^ = 1 (single counting), N^ = 
8,16,24,32 lattices. Linear fits are represented by solid lines, the Tc estimate 
with errors by vertical lines. 

3 String tension and physical scale 

In the previous section by measuring the critical Binder cumulant value we 
obtained high precision critical coupling values for various Nr lattices. In this 
section we measure the Polyakov loop correlator and correspondingly obtain 
the static quark-antiquark potential. The results of the previous section make 
possible to study the scaling of dimcnsionless ratio Td ^fa well into the strong 
coupling regime (/3 = 1.87380). 

To fix the scale we use the critical coupling estimates of the previous sec- 
tion together with the value from the literature [T^ for Nr — 12 lattice. The 
critical temperature T^ = 1/{N^ * a) fixes the lattice spacing in physical units. 
The lattice is then simulated at these Pc values but at different N^- — 3 • N!f 
corresponding to T = Tc/i. At this temperature finite temperature correction 
should be minimal. 

For the string tension measurements we use Liischer-Weisz multilevel algo- 
rithm [ini Uni ■ As was noted in the second reference the one level algorithm is 
computationally preferable. In our simulations the lattice is sliced into layers 
of thickness 2a. On each layer we perform 10 sweeps for iV^/3 = 12 lattice, 100 
sweeps for N^/i = 8 lattices, and 1000 sweeps for A^^/S = 4 and 2 lattices. 
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Figure 8: The Binder cumulant 174 for N-r — 1 (double counting), No- = 
8, 16,24,32 lattices. Linear fits are represented by solid lines, the Tc estimate 
with errors by vertical lines. 



We fit the potential using the following ansatz 



V{f) = Vo-^+af, 
r 



(7) 



where the hats indicate lattice dimensionless observables. 

The fit of the potential is performed in the range f S [2,5 — 9], which ensures 
that no short distance 0(l/r^) or long distance (reduced signal/noise or effects 
from propagation across the periodic boundary condition) artifacts contribute. 
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2.6355(10) 
2.51098(58) 
2.29850(6) 
1.87380(3) 
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0.5505(10) 
0.5826(6) 
0.313(13) 



0.273(35) 
0.2757(12) 
0.3241(7) 
0.256(15) 



0.0153(53) 
0.03232(18) 
0.13312(10) 
0.6285(26) 



0.67(12) 
0.6953(19) 
0.6852(03) 
0.6307(13) 



Table 4: The parameters of the potential fit in lattice units and the string 
tension in (physical) units of Tc for various Nr/'i = N^ lattices. 

In table [4] we collect the data from the potential fits to the ansatz ([7]) . The 
typical goodness of fit Q G [0.42,0.99]. Note that for Nr/i = 2,4,8 lattices 
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Figure 9: The Binder cumulant 34 for N^ — 8, iVg. — 16, 24, 32 lattices. Linear 
fits are represented by solid lines, the T^ estimate with errors by vertical lines. 

we performed the simulations at /3 values from our earlier estimates, which are 
slightly different from the final values presented in the previous section. For the 
errors on the inverse coupling we take either the difference to the final value or 
the final value error, depending on which one is larger. 

In order to compare the string tension a obtained on different lattices, we 
need to convert the lattice observable to physical units. The string tension in 
physical units is ct = a jc? . We can express the lattice spacing in physical units 
through the critical temperature a = ^jiTc ■ ^f)- Therefore we can construct 
a dimensionless observable T^/ \/a = (7V^ • \/5-)~^, which we present in the last 
column of the table. The uncertainties come from the estimates of a and the 
string tension a 

Tr Tr-a , ^ 



Therefore the error of this observable is 




(9) 



where 60 = llAfc/(487r^) = llNc/{2ATr'^) and comes from the scaling of the 
lattice spacing with the coupling in the continuum limii|j. We assume absence 



^Here we consider only the leading term of ^. 
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of correlations between measurements of critical coupling and string tension. 

The data indicates that the string tension values in physical units are con- 
sistent for lattices Nr/?> > 4, although Nr/2> = 12 value obviously needs better 
statistics. The smaller Nt/?, = 2 value is significantly lower than other values. 
This is clearly related to the violation of the hyper-scaling of observables as the 
lattice coupling is increasing. Here we observe the change from the weak cou- 
pling to strong coupling regimes. The scaling window starts around (3 ^ 2.29850 
{Nt/S — 4). It is interesting that the violation of scaling at /3 = 1.87380 is rel- 
atively small (< 10%) for this ratio. Note also that Vq in physical units scales 
like 1/a, therefore Vq is approximately constant ~ 0.5 — 0.6 except for N!f = 2 
where it is significantly lower. 

4 Summary 

We systematically studied N^ = 1,2,4 and 8 finite temperature SU{2) lattice 
gauge theory. The measurement of Polyakov loops in the vicinity of the tran- 
sition point allowed us to study the scaling of the Binder cumulant. We found 
that the critical values of Binder cumulant correspond to the 3D Ising model 
universality class. 



Nr 


f3c 


16 


2.7310(20)* 


12 


2.6355(10)* 


8 


2.5104(2) 


6 


2.4265(30)1' 


4 


2.2991(2) ', 2.2993(3)" 


2 


1.87348(2) 


l(s) 


0.85997(10) 


1(d) 


0.86226(6) 



Table 5: Critical inverse coupling /3c for different N^. lattices: * is from [TH], ^ 
is from [TU], ' is our FSS estimate, while " is obtained from the intersection of 
the largest lattices. 

New high precision estimates of the inverse critical coupling are obtained 
and summarized for various N^- lattices in Tab. [5] together with estimates from 
the literature for lattices where we did not perform measurements. In particular 
we present two formulations for Nt = 1 lattice with results different from the 
previous estimates. 

From the study of static quark-antiquark correlators we extracted the string 
tension and using the critical couplings were able to obtain the dimensionless 
quantity T^j ^fa. This quantity shows small scaling violations for iVr/3 = 2 
lattice (iV^ = 2, /3 = 1.87380) and virtually no violations for larger N^- lattices 
(weaker coupling). 
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